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We study the decay of a two-component Bose-Einstein 
condensate with negative effective interaction energy. With 
a decreasing atom number due to losses, the atom-atom in- 
teraction becomes less important and the system undergoes 
a transition from a bistable Josephson regime to the monos- 
table Rabi regime, displaying oscillations in phase and num- 
ber. We study the equations of motion and derive an analyt- 
ical expression for the oscillation amplitude. A quantum tra- 
jectory simulation reveals that the classical description fails 
for low emission rates, as expected from analytical considera- 
tions. Observation of the proposed effect will provide evidence 
for negative effective interaction. 
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The experimental realization of Bose-Einstein conden- 
sates |lj] offers the possibility of studying quantum matter 
waves. A paradigmatic effect in such systems is the co- 
herent exchange of particles between weakly connected 
states, first predicted for superconductors g. This so- 
called Josephson effect is expected to be displayed also by 
a Bose-Einstein condensate (BEC) in a split trap [^-£|. 
It has been argued that an "internal" Josephson effect 
may also occur in condensates whose atoms can reside in 
two different hyperfine states coupled by Raman transi- 
tions P-pi . Depending on the relative strength of the in- 
teraction between atoms and the hopping between states, 
the system may lie in the Rabi regime, where atoms be- 
have as independent, or in the Josephson regime, where 
interactions are sufficiently strong to create a new col- 
lective behavior B . The coherent exchange of atoms in 
the Rabi regime has been realized between two hyperfine 
states of ^"^Rb atoms (lO|-|l|]. 

It is possible to have two-component condensates in 
which it is energetically favorable for atoms to be in the 
same state. This "ferromagnetic" coupling results from 
an appropriate combination of the three relevant scat- 
tering lengths ||l3| , and does not necessarily require a 
negative scattering length |1J]. On the other hand, atom 
losses by one- and three-body collisions cause the con- 
densate to decay. In this article, we show that the de- 
cay of a two-component condensate with ferromagnetic 
coupling may display spontaneous oscillations due to a 
fundamental dependence of the energy on the relative 
particle number. 

We consider Bose condensed atoms with two internal 
hyperfine levels \A) and \B). The atoms interact via 



AA, BB, and AB elastic collisions. Additionally, a laser 
field induces Raman transitions \A) ^^ \B), generating 
an internal Josephson effect 0,|ll|- If the orbital wave 
functions of states \A) and \B) may be regarded as rigid, 
the system can be adequately described by the two-mode 
Hamiltonian ||6|-p|,[l3|] 
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where the Rabi frequency uiu depends on the detailed 
laser tuning, and where the "charging" energy Ec = 
(uAA + UBB — '2uab)/4: encompasses the various inter- 
actions. The effective interaction constants are given by 
Uij = {AttTi ay/TTT.) J d'^x|(^i(a;)p|(/?j(a;)p, where the mode 
functions '^i{x), i = A,B, are normalized to unity. We 
have assumed uaa = ubb, so that a contribution linear 
in the atom number difference vanishes. In the present 
work we focus on the case where the effective interaction 
energy for the internal Josephson dynamics is negative 
{Ec < 0). This is also the case for a double- well conden- 
sate with a negative scattering length. 

First we analyze the structure of the classical phase 
space. For that in the Heisenberg equations of motion, 
we replace the operators a,b by Y^nje""^-*, ^nse^"^^. 
For the variables z = {ua — nB)/N = 2n/N and = 
{(j^A — 4>b) this yields 
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Hereafter, A = —NEc/2hu!j^, > gives the effective in- 
teraction strength, and the time is measured in units of 
l/cji^;. Equations (|3),(||) can also be obtained from the 
classical nonrigid pendulum Hamiltonian |3l 



H{z,(b,\) = -VT~ 
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where (z, 0) are canonically conjugate coordinates. 

The pendulum Hamiltonian (m obeys the symmetry 
H{z,4>,X) = —H(z,4) + TT,—\) and, thus, a sign change 
of the effective interaction A corresponds to a phase shift 
(/) ^ + TT plus time reversal. This means that the tt- 
states which are found for positive Ec [|l^ at maxima 
of the Hamiltonian, correspond to minima at (/) = for 
negative Ec- 



The classical fixed points of the pendulum Hamiltonian 
can easily be found by linearizing around its extrema. For 
< A < 1 one finds the stable solution 



zo = 0, 



0, r^o = VT^A, 
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where Qq is the frequency of small oscillations around 
the potential minimum. For A > 1 its potential curva- 
ture becomes negative. Thus at A = 1 there is a bifur- 
cation whereby zq becomes unstable and splits into two 
symmetry-related fixed points 



z± = ±v/l-l/A2, 
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As A is proportional to the atom number, a decaying con- 
densate will necessarily cross the bifurcation point A = 1, 
a process which we study below [M . 
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FIG. 1. Classical time evolution of the relative number 
imbalance z{t) (dashed) for 7 = 0.01 and Ao ~ 2. The full 
line denotes the equilibrium value for the corresponding total 
atom number. The inset shows the effective potential (M at 
(p = for the interaction strengths A = 1.5 (dashed), 1 (full 
line), and 0.5 (dotted). 

The bifurcation process can be appreciated in the inset 
of Fig. |l], where H{z, 0, A) is represented for several values 
of A. In the conventional case oi Ec > 0, it has been 
noted [B| that the point where A = EcN/2Tiujji — 1 marks 
the limit between the Rabi (A < 1) and the Josephson 
(A > 1) regimes. Thus we may view the transition from 
(^ to (pi) as the negative Ec version of the crossover from 
an interaction dominated to a hopping dominated regime. 
Once A < 1, the dynamics is controlled by hopping and 
becomes insensitive to the sign of the interaction. 

The master equation that results from the inclusion 
of random losses is such that the classical equations of 
motion (||),(y) remain formally the same, but the effective 
interaction constant acquires a time dependence. 
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due to the exponential decrease of the total atom num- 
ber. We integrate the classical equations of motion 



(g),(||), and (0) with the initial condition z{0) = z+(Ao), 
0(0) = 0, and Aq = 2, i.e., we start from a minimum 
of the pendulum Hamiltonian (Q) in the bistable regime, 
where most of the atoms are in the state \A). Figure |l| 
shows the corresponding time evolution of the number 
imbalance. As the total number decreases due to losses, 
the two minima of the effective potential merge towards a 
single minimum at z = (see inset). The system evolves 
adiabatically until A approaches and crosses the critical 
value of 1. We will see that a departure from adiabatic- 
ity necessarily occurs when A gets sufficiently close to 1. 
After crossing the bistability threshold the system finds 
itself distant from the new equilibrium point and begins 
to oscillate spontaneously around it. These oscillations 
initiate in the crossover region between the Josephson 
and the Rabi regimes. 

The oscillation amplitude can be estimated analyti- 
cally if one assumes (i) that outside the transition region 
A « 1 the effective potential changes sufficiently slowly 
and (ii) that the crossover region is traversed fast enough 
for the process to be described as a sudden change in the 
potential. 

For A > 1 the condition for adiabaticity reads 
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while later, for A < 1, the condition is 
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For A « 1 the r.h.s. of the inequalities (^) and (Q) become 
7/[2(A - l)]3/2 and 7/2(1 - A)^/^, respectively. This in- 
dicates that, sufficiently close to the bifurcation point 
A = 1, adiabaticity breaks down. 

Initially, the system follows the potential minimum 
z+(A(t)) adiabatically until the inequality (||) is violated 
for the ffi'st time. This occurs for 



z = z+(A) 
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where Awl has been assumed. As A crosses the transi- 
tion region, we may assume that there is a sudden change 
of the effective potential in the sense that z and (p are ef- 
fectively constant. Then the adiabatic condition, now 
reading as in (0), is recovered for A sufficiently below 1. 
The resulting dynamics can be conveniently studied 
by noting that, for a harmonic oscillator, the action J = 
2ttE/V.o is an adiabatic invariant p7l]. T he energy E = 
\{(jp -f rjpZ^) is estimated from Eq. (|10D leading to J = 
27r51oz^ — 7r7. Thus the instantaneous amplitude reads 
[7/i7o(''')]^^^ and the imbalance oscillates like 



z{t) w v^7/rio(T)cos[rio(T)r], 



(11) 



where r is the time elapsed after the oscillation sets in, 
and ^^{t) = (1 — e"''''^)^/^ is the instantaneous oscilla- 
tion frequency. For sufficiently small 7, ^o{t) — (7t)^/2 



which yields, for the time r^ at which the fc-th max- 
imum of the oscillation happens, the condition 2Trk = 

^o(T'fc)Tfc ~ 1^^"^ '''k ■ "^^^ corresponding amplitudes 
therefore read 



Zk « (7V2^fc)'/'- 



(12) 



We have numerically checked the values of Zk for 1 < fc < 
5 in the range 10~^ < 7 < 10~^, and have found good 
agreement with this analytical prediction. 

The classical treatment above requires a smooth po- 
tential, such that its curvature is basically constant on 
the scale of the width Az of the wave packet. While 
approaching the transition point A = 1, the potential 
becomes flatter and, thus, the stationary wave packet 
spreads. At the same time the distance |z+ — z_| be- 
tween the two potential minima decreases. As long as 



Az < \z+ 



(13) 



the classical description should be adequate. When in- 
equality (13) is violated, a quantum description is nec- 
essary. For that we take the variables z and 4> to obey 
the commutation relation [n, e*"^] = ze*"^, which for small 
phase uncertainty amounts to [z, 0] = 2i/N. For the 
Hamiltonian M) in the harmonic approximation, this 
yields the number r.m.s. (Az)^ = l/Nfl^. 

When A ^ 1+, the adiabaticity condition (m reads 
[2(A - l)]3/2 > 7, while (|l|) amounts to [2(A - 1)]^/^ > 
l/AN. We expect the classical description of the passage 
through the bifurcation to be correct if, as A decreases in 
the bistable region, adiabaticity breaks down before clas- 
sicality does 18 1 . This is the case if A^7 > 1, i.e., if losses 
are sufficiently intense that more than one atom is lost 
in a typical Rabi cycle. On the contrary, we expect the 
classical predictions to fail for iV7 < 1. A similar conclu- 
sion is reached for the regime A < 1 if we replace ( |l3[ ) by 
the condition that the width {NQq)'-^^^ of the ground 
state wave function is smaller than the scale y/l — A be- 
low which the effective potential behaves as harmonic. 
Below we see that these expectations are confirmed by 
the numerical simulation. 

To analyze the quantum dynamics of a single experi- 
ment, we employ the quantum jump method p9|| where 
the time evolution over a sufficiently small time interval 
At is simulated as follows. If an atom is lost from conden- 
sate A (B), the wave function changes as \ip{t + At)) — > 
Q\ip{t)), Q = a,b. These losses occur with the proba- 
bilities pqit.At) = ^At{i}]{t)\Q'^Q\i}}{t)). If no loss takes 
place, the system propagates form t to t + At with the 
effective Hamiltonian i?eff — H — ihrf{a^ a ^-h^h) / 2. Since 
the propagation with i?cff as well as the atom losses are 
non-unitary, the wave function has to be normalized after 
each time step. Note that in the case of two identical loss 
rates the total atom number commutes also with the ef- 
fective Hamiltonian and terms in iJcff which depend only 
on the total atom number can still be neglected. 



We start our simulation in the bistable region with 
Ao = 2. Thus the initial state is the "local ground state" 
at the potential minimum z_|_ and satisfies condition (^3). 
Here we make the very reasonable assumption that the 
parity symmetry is effectively broken due to the large 
value of iVo, i.e., that any Schrodinger cat state has 
evolved into a mixture due to tiny perturbations from 
the environment. 

Figure 0(a) depicts the quantum mechanical time evo- 
lution for the same parameters as used in Fig. n^. Both 
plots agree quite well, which indicates that the system be- 
haves classically. However, for a smaller value of the de- 
cay rate [Fig. |^(b)] we find substantial differences: (i) the 
oscillation starts well before the transition to the monos- 
table regime and (ii) its amplitude is larger than in the 
corresponding classical case (see Fig. ^. This is so be- 
cause the finite width of the wave packet allows to cross 
the point z = when there is still a finite barrier. 
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FIG. 2. Quantum mechanical time evolution for 

iVo = 1000 and Ao = 2. The decay rate is 7 = 10"^ (a), 
10~* (b). The steps in the stationary value z+ are due to the 
discreteness of the total atom number. 

A comparison of the classical and the quantum me- 
chanical oscillation amplitude (Fig. g) shows that the 
quantum dynamics agrees to the classical dynamics for 
sufficiently large values of the decay rate 7. We find 
that the quantum mechanical oscillation amplitude is 
larger than its classical counterpart for 7 < 3 x 10^^^, 



i.e., roughly at the value for which the classical descrip- 
tion was predicted to fail. Then the amplitude is com- 
parable to the phase uncertainty and therefore difficult 
to observe. A decay rate Nqj « 10 seems optimal for 
the experimental observation, since then there are still 
sufficiently many oscillations before the condensate has 
decayed and the amplitude is still larger than the uncer- 
tainty. 
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FIG. 3. Amplitude of the quantum mechanical oscillations 
(circles) and its quantum uncertainty (boxes) for A'o = 1000 
and Ao = 2. Each data point results from an average over 20 
simulation runs. The full line displays the classical result (113) 
for k — 1. 

Let us assume that the rigid wave functions ipa,b{x) 
are given by the Thomas-Fermi (TF) approximation for 
1000 ^^Rb atoms in an isotropic trap with frequency 
100 Hz. This results in a TF radius of 4/im and a 
harmonic oscillator length of 1 fim. An effective neg- 
ative scattering length Ocfr = uaa + o,bb — '^ciab ~ 
-qaa/WO « -0.5 nm yields Ec/h « -0.3 s"^ Then, 
for A^o = 10'', to have Ao = 2 requires a Rabi frequency 
ujR « 27r • 10 Hz, the "optimum value" for the decay rate 
corresponds to a condensate life time of {•^ujj^)^^ ^ 1.5 s. 

In summary, we have shown that a two-component 
condensate with negative interaction energy (Ec < 0) 
will exhibit spontaneous oscillations in phase and num- 
ber during its decay due to atom losses. These oscil- 
lations have their root in the fundamental dependence 
of the Josephson coupling energy on the relative boson 
number. As the atom number decreases, interaction be- 
tween atoms becomes less important than hopping be- 
tween states. Thus the system evolves from a Josephson 
regime, which is bistable due to the negative sign of Ec, to 
a monostable Rabi regime dominated by hopping. A clas- 
sical description of the oscillatory decay dynamics gives 
a good quantitative account of the decay process for high 
enough loss rates. The observation of this effect will pro- 
vide unambiguous evidence of the Josephson effect with 
negative Ec, by identifying the existence of a bifurcation 
that is manifested through the spontaneous appearance 



of Josephson-Rabi oscillations. By accurately monitor- 
ing the decay process, it will be possible to quantitatively 
check our current understanding of the internal Joseph- 
son dynamics of two-component condensates. 
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